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Abstract 

In contrast to one-component Bose-Einstein condensate case, the vortices in two-component 
condensate can have various complicated structures. The vortices in a space-homogeneous 
Bose-Einstein condensate have been studied in this paper. It is shown that the vortex struc¬ 
ture is described by three dimensionless parameters. This is totally different from the usual 
one-component condensate case, where an isolated vortex is described by a parameter-less 
dimensionless equation. The two-conrponent vortex structure strongly depends on the sign of 
“interaction” constant of the components. A few types of vortices with different qualitative 
structure are explored. We show that the super-density vortices can exist, when the “interac¬ 
tion” constant is positive. The super-density vortices have the near-axis density greater than 
the equilibrium density of a homogeneous space Bose-Einstein condensate. We also show that 
the vortices with opposite direction of the condensate component rotation near the axis and 
far off the axis can exist. 

PACS: 67.85.Fg 

Key words: Bose Einstein multi-component condensate, Gross-Pitaevskii equations, quantum 
vortex. 

In [1] [2] the quantized vortex filament in a weakly nonideal Bose gas was investigated. The 
structure of the vortex is on the axis of the vortex the following density of a Bose-Einstein con¬ 
densate (BEC) vanishes and increases monotonically to the equilibrium value of the density of 
a spatially homogeneous system at infinity. The characteristic distance over which the density 
varies near the vortex axis is determined by the correlation length, which depends on the mass of 
the particle, the equilibrium density and interparticle interaction constant. The correlation length 
determines the size of the vortex core. The equation describing the vortex in a single-component 
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system, if it is written in dimensionless form, is universal and does not contain any parameters, 
characterizing the system [3j. Consequently, the structure of the vortices in all single-component 
systems is the same. The influence of nonlocal effects on the vortex structure was seen in [3j. 

Creating a BEC in atomic gases of alkaline elements contained in magnetic or laser traps mum 
IB], led to experimental study of the vortices in spatially inhomogeneous conditions [9,TQJ UH! .12, 13] . 
Were examined not only single-component systems, but has been implemented the ability to create 
condensates mixture of two or more types of particles and vortex states in these more complex 
systems [in 112113]. Although the two-component BEC in trap with vortices have been studied 
theoretically [2 [8] [H, 15] HHi H2, IB] [19], a detailed analysis of the structure of the vortex in 
such systems, depending on the characteristics of the environment, as we know, is absent for the 
spatially homogeneous case. 

In this paper we study the structure of vortices in the spatially homogeneous Bose-Einstein con¬ 
densate of two kinds of particles with different masses. It is shown that in this two-component case 
a structure of the vortex depends on the three-dimensional parameters characterizing the system. 
Changes in the density of each component, the total density and mass flux density as a function 
of distance from the axis of the vortex is determined by the values of these parameters. Attention 
is drawn that the considered system differs significantly from the one-component case where the 
vortex is described by a universal nonlinear equation in dimensionless form which does not contain 
the characteristics of the system mmm. Some possible qualitatively different vortex structure 
with different sets of parameters characterizing the two-component condensate is explored. It is 
shown that the structure of the vortex depends on the sign of “interaction” component constant. 
In particular, the existence of a sort of “super-density vortices” is possible in case with positive 
value of the coupling constant. In such vortices total density near the axis may exceed the equilib¬ 
rium density spatially homogeneous condensate. This is impossible in single-component case. It 
is also shown that there are vortices in which the condensate near the axis and at large distances 
from the axis of the vortex is rotated in opposite directions. 

In the Erst section of the paper the system of the Gross-Pitaevskii (GP) equations for a two- 
component BEC is presented in dimensionless symmetrical form. The system contains five dimen¬ 
sionless parameters, among them only three parameters are independent. In the second section 
the GP equations for macroscopic wave functions absolute values are written in the cylindrical 
coordinate system in the presence of vortex motion with two quantized circulation. In the third 
section we investigate the asymptotics of solutions at large distances from the vortex axis. It is 
shown that the solution of the equations may tend to its asymptotic value from above and below, 
thus under certain conditions the total number density near the axis of the vortex can be greater 
than the equilibrium density at infinity. Some qualitatively different vortex structures are numer¬ 
ically studied in the fourth section. We consider the cases when the circulation is different from 
zero in one of the components, or in both of them. In the fifth section we analyze the distribution 
of the mass flow density in the vortex filaments. 
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1 Gross-Pitaevskii equation for two-component system 

Bose Einstein condensate, consisting of two particle species with mass mi, m 2 and densities ni, n 2 
is described by stationary Gross-Pitaevskii equations for macroscopic wave functions 


- 77 — A$i ( r ) + 9n®i (r) |<hi (r)| 2 + 9 i 2 $i (r) |<f> 2 0)| 2 

Zuij\ 

h 2 

-- —A<h 2 (r) + 922^2 (r) |$2 (r)| 2 + 92i$2 (r) \®i (r)| 2 
Zrri2 

The particle number for each species is defined by relation 


hi$i (r), 


/i 2$2 (r). 


( 1 ) 

( 2 ) 


N l = J \^(r)\ 2 dr, (i — 1,2). (3) 

We use the index i for denoting particle species and condensate components. Complex functions 
can be represented through the absolute values pi and phases Xii writing them in the form of 
<f>j (r) = pi (r) exp [i\i (r)]. Then the particle number density rii (r) and the particle velosity v* (r) 
are defined by the absolute value and the phase of complex function 


0) = P 2 i ( r ), 


v * (r) 


h 

rrii 


Vxi (r) • 


In homogeneous state the chemical potentials pi, p 2 are 


(4) 


hi — 9ll n 10 + 9l2’ n 20, p2 ~ 922 n 20 + .921 10- (5) 

We assume that the stability condition for homogeneously system is satisfied mi 


.9 ii > 0 , 922 > 0 , 9n922 — .9 12 > 0 . 

We use the dimensionless functions for convenience: i])\ = $ 1 /ybi 10 , V ’2 = $ 2 / y/n 2 Q. 
The healing lengths for the first and the second condensates are as follows 


a = 


h 

\/2mi9iPTo’ 


6 = 


h 

\/2m 2 g 2 2n 20 


( 6 ) 


(7) 


Dehne effective healing leangth £ = . Then the dimensionless coordinates will be r = r/( 

and A = £ -2 A. Hereafter the tilde symbol will be omitted. Dimensionless form of equations (JT|) 
and ([2]) reads 

-AVh (r) + Aid^i (r)| 2 - l)' 0 i (r) + uB 1 (^2 (r)| 2 - l)^i (r) = 0, . , 

-Ai/i 2 (r) + H 2 (|^2 (r)| 2 - 1)^2 (r) + (r)| 2 - l)-02 (r) = 0. 

The system of equations (jSJ) has a symmetric form. It has five dimensionless coefficients: Ai, A 2 , Bi, 
B 2 and u = 912/ \j9u922- Coefficients A t , Bi are not independent. They are expressed by two 
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dimensionless positive parameters, one of them is defined by ratio of healing lengths ( 2 = (C, 2 /fi) 2 , 
and the second of them is defined by ratio of densities and interaction constants in each component 
7] = ( 7120 / 7110 ) \J 922 ) 9 ii- Thus the solution of the system of equations (JT|) and ([2]) depend on the 
following three independent dimensionless parameters 

>2 \ 2 ri 20 1922 912 

C = ( 7“ ) > V = — \ —, u = (9) 

Vsi / n w V 9i 1 \J9w922 

The first two parameters in ([2]) are positive but the third parameter satisfies inequality 7t 2 < 1. 
The coefficients in ([S]) are expressed through parameters (J2J) : 


^i = l + C 2 > ^2 = 1 + C 2 > 

B\ = Air), B 2 = A 2V - 1 . { 1 

In single component case the dimensionless GP equation has not the parameters of medium 
and thus it is the universal equation. In two component condensate the coefficients of system are 
defined by three dimensionless independent parameters (TUT) . The vortex has different structure 
depending on the values of this coefficients. 


2 Gross-Pitaevskii equation for vortex 


We seek for solution of the set of equations ([HD with the symmetry about the axis of 2 in the form 
(r, Lp ) = fi (r) exp [i\i (</?)], where r, (p are the cilindrical coordinates. 

In this case the velocities of the components have the form 


v* = 


h dXi (p) 

777 , iT dip 


where is the unit vector tangent to the circle with center on the z-axis. 
The set of equations for functions /* (r) from ([HD are: 


(ii) 


( 12 ) 

(13) 

We will study the solutions of this equations with the following boundary conditions: 
if l\ > 1 we take fi (0) = 0 and fi ( 00 ) = 1, when li = 0 we take dfi/dr( 0) = 0 and fi(oo) = 1. 

The total density of the particle number in the vortex is 


d 2 fi 
dr 2 

fh 

dr 2 


1 dfi l\ . 2 

+ ^ /i+ i(/i 

+ 1 /2 + am 


l)/2 + uB 2 (fi — l)f 2 


0, 

0. 


n (r) = n w fl (r) + n 2Q f 2 (r). (14) 

The total particle number density is equal to the equilibrium density away from the axis of the 
vortex n ( 00 ) = n 0 = nio + n 2 o- 
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We take the phase of macroscopic wavefunction in the form (t 5 ) — hp, where U = 0, ±1, ±2,.... 
The sign of U determines the direction of vortex rotation. In a two-component system the vortex 
motion in each component has its own quantized circulation: 


r, = 


■dl = 


2i\hL 


rru 


(15) 


An interesting feature of a two-component condensate is the possibility of vortex existence in 
one component, while the other component remains irrotational with the circulation equal to zero. 
In this case l\ = ±1 and Z 2 — 0, or li — 0 and Z 2 = ±1. 

Another significant feature is manifested in the rotatable component in opposite directions. In 
particular, l\ = +1 and I 2 = —1 corresponds to this case. Mass flow density in the vortex filament 
of two-component condensate is defined by the relation 


J = “ [^ 10/1 ( r ) + kri 2 ofi (r)] e v = J 4> ( r ) e v ( 16 ) 

With this in mind, the angular momentum of the vortex L = f [r x J] dr in a two-component 
BEC is given by 


L = h{l 1 N 1 + l 2 N 2 )e g , (17) 

where e z is a unit vector along the axis of the vortex, and N t = 2nnioL v fo° ff ( r ) r ^ r i s the total 
density of particles of species i, L v is the length of the vortex. Note that the density distribution 
of the component does not depend on the sign of circulation, but the angular momentum (fTTlh of 
course, depends on it. If the rotation of component is the opposite, and the condition Zinio+Z 2^20 = 
0 is satisfied, the total angular momentum of the vortex vanishes. 

3 Asymtotic solutions of equations for vortex 

We consider the asymptotic behavior of solutions of (!T2]) . (lT3j) . The solutions of this equations have 
the following form at large distances r>l: /j « 1 + aj/r 2 where 

{l\M ~ urjl^Ax) {I2.A1 - ur]^ l llA 2 ) n oN 

Ul ~ 2 AiA 2 (1 — u 2 ) a2 ~ 2A\A 2 (1 — u 2 ) ' ( 

The denominator in these formulas is always positive and the numerator can have different 
signs. Depending on the sign of the parameter a, : ,the absolute value of macroscopic wave function 
may approach to its equilibrium value (unity) at infinity from above or from below. Obviously, 
if u < 0, then both parameters a* are negative. If u > 0, then it is possible that one of these 
parameters is positive, and the other parameter is negative. Suppose l\ > 0, l\ > 0. Then, for 
u > 0 the following cases can be implemented: 

1) if A 2 /Ai <ur] then cq >0, a 2 < 0, 

2) if urj (l\fl 1 ) < A 2 /A 1 Ku^r] (l\/l f), then a\ < 0, a 2 < 0, 

3) if A 2 /A 1 >u^ 1 r] {l\fl 1 ), then a\ <0, 02 > 0. 
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If circulation is absent in one of the components, for example, in the second component 
l 2 = 0, then ai = —l\A 2 /2A 1 A 2 (1 — w 2 ), a 2 = u , g~ l l\A 2 /2A 1 A 2 (1 — u 2 ). 

In this case a\ < 0 and the sign of a 2 can be either positive or negative, and coincides with the 
sign of the constant of the interaction between the components u = g\ 2 /yjgug 22 . 

According to (fT4|) and (TT8|i the behavior of the total density at large distances defined by the 
formula 


n (r) _ 2a D 

no 7i 0 r 2 

where 


(19) 


&D — Oinio + fl-2^20 


We introduce the notation 


[(ni 0 / 2 A 2 + n 2 ol\Ai) - u (ni 0 't]l 2 Ai + n 2 pg l l\A 2 )\ 
2A,A 2 (1 - u 2 ) 


( 20 ) 


(nip/f A 2 + n 2 pl 2 Ai) 

{n 2 or]-HlA 2 + n w ril%Ai)' 

When u = u* the parameter old (120]) is equal to zero. In case of positive parameter w* > 1, for any 
admissible values of u we have ao < 0 and, therefore, the total density tends to its equilibrium 
value at infinity from below. An interesting possibility arises if w* < 1. It is possible to have 
u > u* and old > 0 in this case. This means that the total density approaches its equilibrium value 
at infinity from above, i.e. density near the vortex axis is greater than the equilibrium density at 
infinity. We emphasize that the possibility of the existence of such “vortices with super-density” 
near the axis is a feature of the multicomponent condensation. When the circulation in the second 
component is absent (l 2 = 0), an has the simple form 


_ l{ (n w - ut) 1 n 20 ) 
a ° ~ 2A 1 (1 - u 2 ) 


( 22 ) 


Coefficient (122]) is positive if the condition n w > ur]~ l n 2 o is satisfied. This condition is equivalent 
to g 22 > g\ 2 in accordance with ([9]). In this case, the total density near the vortex axis is greater 
than the equilibrium value of the density at infinity. 


4 The vortex structure on the different signs of the inter¬ 
component interaction constant 

We investigate numerically the structure of the vortex at different sign and magnitude of the 
parameter u, which determines the “interaction” of BEC component. Consider first the case when 
the vortex motion in the second component is missing, so the l 2 = 0. Put in numerical calculations 
that C — V — 1, thus, A\ = A 2 = B\ = B 2 = 2, and w* = n\o/n 2 o (121]) . Inter-component 
interaction parameter can vary in the range — 1 < u < 1. 
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Numerical solutions of the equations (|12p . (I13p and the behavior of the total density as a 
function of distance from the vortex axis at u — —0.7, u = —0.98, u = 0.9 and u = 0.6666 are 
shown in Figures Q] and [2j In the plots we follow the notation with dimensionless r. The case 
u = —0.7 is shown in Figure QJi. The macroscopic wave-function absolute value for component 
with circulation vanishes on the axis. In component, where there is no circulation, it takes a finite 
value less than unity. The wave functions of the two components and the relative total density 
increases monotonically with distance, approaching the equilibrium values. The total density on 
the vortex axis, in contrast to the usual single-component BEC vortex, does not vanish. 

Figure lb shows a similar dependence at u — —0.98. In this case, we are near the stability 
limit of the system (see ([6])). We see that the qualitative behavior of the total density is the same 
as that of the single-component condensate, but the value of the total density on the axis of the 
vortex is nonzero. 

Figure 2a shows graphs at u = 0.9, w* = 2/3 and u* = 1. In the case u* = 2/3, in contrast 
to the previous case, the magnitude of the wave function components without circulation on the 
vortex axis is greater than its value at infinity, and it decreases monotonically with increasing 
distance, approaching unity from above. Absolute value of the wave function components with 
vortex motion increases monotonically with distance from zero to the equilibrium value. 

The behavior of the total density in this case may be different: if u > w* it decreases mono¬ 
tonically from a value greater than the equilibrium value (upper solid curve in Fig. [2]r), and at 
u < w* density increases monotonically from the smaller unit on the axis (the lower solid curve in 
Figure [2]r). 

Some interesting features in the structure of vortices in two-component condensate are found 
for certain relations between the parameters u and w*. The case u = 0.6666, u* = 0.5625 is shown 
in Figure [2b- Here the total density remains almost unperturbed, though the circulation is present 
in the system. 

The numerically obtained vortex solutions with finite density on the axis and vortex solution in 
which the density on the axis greater than the equilibrium are qualitatively different from the usual 
vortices in one-component BEC. Numerical calculations for a smaller absolute value of constant 
interaction |u| show that the qualitative behavior of solutions of the equation GP and the total 
density does not change in this case. When the interaction between the components is reduced, the 
distribution of the wave function in the irrotational component becomes more uniform, approaching 
the uniform distribution of / 2 = 1. 

Consider the case, where two circulations are associated with a vortex, l\ — — 1. Here, 

as above, we believe that £ = 1, parameter r] can have any positive value. Thus, A\ = A 2 = 2 
B\ = 2 r], B 2 = 2 r]~ l . Inter-component interaction parameter can vary in the range —1 < u < 1. 
Here, as in the above case, there may be situations when the density of the components approach 
their equilibrium values at infinity from below and from above, depending on the sign of constant 
^ (ITS]) . The following cases are possible: 

1) if u > 7] we have oq < 0, a 2 > 0; 

2) if i] > 1/iiwe have cq > 0, a 2 < 0; 

3) if u < r] < 1/u we have a± <0, a 2 < 0. 

Figure [3] shows the dependence of the solutions of the GP equations at the parameters, when one 
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Figure 1: The solutions of GP equations for f\ (r), / 2 (r) and for total density n(r ) in the case 
l\ = 1, Z 2 = 0, w* = nio/n -20 = 2/3: a) w = —0.7; b) u = —0.98. 


of the solutions tend to the value at infinity from below, and the other solution tends from above. 
If the parameter u is less than the critical value of w* (OTj) (Fig. [3k), the total density increases 
monotonically with increasing distance from the axis. 

If u > w*, then the vortex structure is realized in a qualitatively different manner from the 
structure of the vortex in the one-component case (Fig. [Hb). The density increases from zero on 
the axis to a maximum value that exceeds the density at infinity and then monotonically decreases 
with the distance to equilibrium density. In this case vortex has the form of tube nearly empty 
within and limited by solid wall. 

Thus, it is shown that there may be two types of “super-density vortices”. In one case, the 
density on the axis of the vortex exceeds the equilibrium density and decreases monotonically with 
distance (the upper solid curve in Figure [2k). In another case, having a density equal to zero at 
the axis, the density reaches the maximum value, greater than the equilibrium density, at a certain 
distance of the order of the correlation length (Fig. [3k>). 


5 Vortex mass flux density 


It is interesting to consider how the mass flux density varies with distance from the vortex axis. 
The mass flux density is defined by the formula ffIBD . which is conveniently written in dimensionless 
form 












(b) 


Figure 2: The solutions of GP equations for f\ (r), f 2 (r) and for total density n(r ) in the case 
l\ = l,/ 2 = 0: a) w = 0.9; b) u = 0.6666. The upper solid curve in Figure (a) shows the behavior 
of the total density n if w* = nio/n 2 o = 2/3 < 0.9. The lower curve corresponds to the parameters 
■u* = nio/n 2 o = 1 > 0.9. The Figure (b) shows the same function for parameters: u = 0.6666, 
w* = 0.5625. 
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(a) 


(b) 


Figure 3: The solutions of GP equations for /i (r), / 2 (r) and for the total density n (r) if 
l\ — l\ — 1: a) u — 0.7, b) u = 0.9 (£ = 1, r/ = 2, u* = 0.8, nio = n 2 o). 


^( r ) 

J<po 


, n 10 ,2 / \ , 7 n 20 ,. 2 

*i —fi (r) + h — / 2 

n 0 n 0 


(23) 


where J^ 0 = 7m 0 /^. In the case r€l ^ 0, ^ 2 7 ^ 0 the solution of equations (fT2]h (fl3l) vanish 
under the law /i ~ Cir, / 2 ~ C 2 r. Finding the expression for the mass flux density at small and 
large distances using the condition that r > 1 we have j\ ~ / 2 ~ 1: 




7 0 r, (r < 1); 


^ ( r ) 7oo , 

~ —, (r > 1). 

r 


(24) 


• 700 < 7^0 

Direction of rotation near the vortex axis and at large distances are determined by the signs of 
the coefficients 


mo ^ 2 , , ri2o C 2 


7 o = h — C 1 + Z 2 — 
n 0 n 0 


2t 


Too 


— ; 77,10 I ; 77,20 

= h -r h -• 

n 0 no 


(25) 


The coefficients Cf, C 2 can not be calculated analytically, but the numerical calculations show 
that the coefficients (l25|i may have different signs, when the signs of U and Z 2 are opposite. This 
means that in this case the condensate in the vortex near the axis and at large distances is rotated 
in opposite directions. Figure Q] shows the dependence of the mass flow density in the vortex of 
the distance from its axis. The solid curve corresponds to the case when the signs l\ and Z 2 are 
the same (l\ = — In this case, the distribution of the mass flow is not qualitatively different 

from the distribution in one-component BEC (dotted curve). 
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(a) (b) 

Figure 4: Full mass flux density as a function of distance from the vortex axis. The ordinate value 
is (r) / J(p() , where J^ 0 = hn 0 /£. The solid curves correspond to the same sign b and 1-2, and the 
dashed curves correspond to the different sign. Figure (a) shows a case of the same direction of 
vortex rotation. Figure (b) shows the case of the direction change of vortex rotation. The dotted 
curve refers to the usual one-component case. 
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The most interesting case is, when the signs l\ and Z 2 are opposite. Dependence of the mass 
flow density in the vortex on the distance from the axis is shown by dashed lines in this case. 
The dependence on Figure Hk is calculated for a set of parameters: n 10 /n 2 o = 3/2, £ = \J 3/2, 
7 ] = 4/3, u = 0.7. The rotation of the condensate at any distance from the axis occurs in the 
same direction in this case. When parameter set is nio/n 2 o = 2/3, £ = y/2/3, r/ = 3, u = 0.7 
the situation with the condensate rotation near the vortex axis and at large distances in opposite 
directions is realized (Figure Hb). 

Note that some of the effects of quantized vortices with opposite rotation in the two-component 
condensates in atomic gases have also been studied in a recent paper m- 

6 Conclusion 

In the one-component BEC the vortex structure is universal and does not depend on the char¬ 
acteristics of the medium. There is a kind of a law of corresponding states, consisting in the 
fact that the equation for vortex, written in dimensionless form does not contain any parameters 
characterizing the system. This does not hold for the condensates, consisting from particles of 
different species. In particular, a two-component condensate vortices are characterized by three 
independent dimensionless parameters and, depending on the ratio between them, may have a 
different structure. 

Some possible vortex structure bearing one and two quantized circulations in the two component 
BEC are discussed in this paper. Numerical solutions of GP equations are shown at different 
parameters. It is shown that the density distribution in the vortex depends on the sign of the 
parameter describing the interaction between the components. For positive parameter, there are 
solutions with density near the vortex axis greater than the equilibrium density at infinity (“super¬ 
density vortices”). If the vortex has one quantum circulation, its density on the axis can reach a 
maximum value exceeding the equilibrium density. For vortices with two quantized circulations 
the density on the axis is zero and the maximum density that exceeds the equilibrium value is 
achieved at a distance from the axis of the correlation length order. 

Another interesting feature of the vortices discussed in this paper is the ability to change the 
direction of rotation within the vortex. It is shown that vortices with opposite direction of rotation 
near the vortex axis and away from it can exist. The asymptotic behavior of vortex solutions over 
long distances is studied. It is shown that their analysis allows qualitative conclusions about the 
structure of the vortex. 
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